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ABSTRACT 

We study the density structure of the candidate contracting starless core 
L694-2 using 1.3 mm dust continuum observations from the IRAM Plateau de 
Bure Interferometer and the Berkeley-IUinois-Maryland Array, which probe spa- 
tial scales from 10000 AU to 500 AU. The long baseline PdBI observations detect 
no emission from the core, and limit the maximum contamination from a com- 
pact component < 2.7 mJy. The flux limit corresponds to a very small disk 
mass, Mdisk ^ 5 x 10"^ M© (60 K/Tdisk), and bolsters the "starless" interpre- 
tation of the L694-2 core. The shorter baseline BIMA data are compared to a 
series of density models using a physically motivated temperature distribution 
with a central minimum. This analysis provides clear evidence for a turn-over 
from the steep density profile observed in the outer regions in dust extinction to 
substantially more shallow behavior in the inner regions (< 7500 AU). The best 
fit Bonnor-Ebert, Plummer-like, broken power law, and end-on cylinder models 
produce very similar flattened profiles and cannot be distinguished. We quantify 
the sensitivity of the inferred structure to various uncertainties, including the 
temperature distribution, the accuracy of the central position, and the presence 
of a weak unresolved central component. The largest uncertainty comes from the 
temperature assumption; an isothermal model modifies the best fit parameters 
by ~ 2(7, with the inferred density profiles more shallow. Dust emission and ex- 
tinction profiles are reproduced by an embedded isothermal cylinder with scale 
height H = 13.5" inclined at a small angle to the line of sight. The turn-over 
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observed in the L694-2 density distribution suggests that pressure forces still sup- 
port the core, and that it has not fully relaxed as in the inside-out collapse model, 
despite the extended inward motions inferred from molecular line observations 
(Lee, Myers, & Tafalla 2001). In the context of the cylindrical density model, 
these inward motions may represent the contraction of a prolate core along its 
major axis. 

Subject headings: ISM: globules — ISM: individual (L694) — radio continuum: 
ISM — stars: formation 



1. Introduction 

The observed properties of dense cores (e.g. Benson & Myers 1989) form the basis of the 
standard model of isolated star formation. In this model, the "starless" dense core represents 
the earhest identifiable stage of the star formation process. The physical conditions in 
this early stage have a profound impact on the evolution of protostars towards the main 
sequence. The initial density structure, particularly in the innermost regions, affects the 
collapse dynamics and the time dependence of the mass accretion rate and therefore many 
of the observable properties of protostars, including luminosity. 

A quantitative understanding of the collapse dynamics has been hindered by uncertain 
knowledge of appropriate initial conditions (Andre, Ward-Thompson & Barsony 2000). In 
the popular theory of "inside-out" collapse, a spherical starless core loses turbulent and 
magnetic support and relaxes to a balance between gravity and thermal pressure, an 
density distribution is estabhshed and the core collapses from the inside-out with a constant 
mass accretion rate (Shu 1977). However, if collapse begins before the density distribution 
fully relaxes, then a central region of relatively constant density remains and the mass 
accretion rate is an order of magnitude larger at early times (Foster & Chevalier 1993). 
This phenomenon has been identified with the youngest "Class 0" protostars, which exhibit 
especially powerful outffows (Henriksen, Andre & Bontemps 1997; Andre, Ward-Thompson 
& Barsony 2000). Better observations of starless cores are needed to determine the initial 
conditions. 



^ Based on observations carried out with the IRAM Plateau do Bure Interferometer and the Berkeley 
Illinois Maryland Array. IRAM is supported by INSU/CNRS (France), MPG (Germany) and ION (Spain). 
The BIMA array is operated by the Berkeley- Illinois-Maryland Association under funding from the National 
Science Foundation. 
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The L694-2 dense core is one of several "strong" infall candidates among starless cores, 
based on observations of molecular line profiles with redshiftcd self-absorption in a systematic 
search of more than 200 targets (Lee, Myers & Tafalla 1999). Molecular line maps of these 
objects provide strong evidence of inward motions, with speed ~ 0.1 km s~^ over a radius of 
~ 0.1 pc (Lcc, Myers & Tafalla 2001). The physical basis for these motions is unclear. The 
speeds are subsonic and may be associated with condensation through ambipolar diffusion, 
or perhaps a magnetically diluted gravitational collapse (Ciolek & Basu 2000). Alternatively, 
pressure driven motions due to the dissipation of turbulence may be responsible (Myers & 
Lazarian 1998). 

Observations of dust column density provide a means to infer the structure and dynam- 
ical state of a dense core. Harvey et al. (2003b) describe near-infrared extinction measure- 
ments toward L694-2 that show a very steep density profile in the region r = 30" to 83" 
(0.036 to 0.1 pc, or 7500 to 20000 AU), where background stars were detectable for extinction 
measurements. The extinction data can be fit by a number of models that produce nearly 
degenerate column density profiles, including a simple power law with index p = 2.6 ± 0.2, a 
supercritical Bonnor-Ebert sphere with dimensionless outer radius .^max = 25 ±3, or a nearly 
end-on isothermal cylinder with scale height H = 13.5" ± 1.5". The latter two models pro- 
vide a physical description of the core that represents an unstable configuration of material 
that is consistent with the infalling motions inferred from molecular spectral line profiles. 
In particular, the slightly tilted cylindrical model that provides a basis for interpreting the 
asymmetry of the L694-2 core should be highly unstable to gravitational collapse along its 
axis (close to the line of sight) , consistent with the strength and orientation of the observed 
velocity structure. The model also provides a framework for understanding the very high 
extinctions observed in L694-2, namely that the core has prolate structure with the full ex- 
tent (mass) of the core being missed by analysis that assumes symmetry between the profile 
of the core in the plane of the sky, and the profile of the core along the line-of-sight. 

Observations of long wavelength dust emission provide another means of probing density 
structure. This technique is almost as direct as observations of dust extinction. The intensity 
of the optically thin emission provides an integral along the line-of-sight of the product of the 
density, temperature, and opacity of the dust. The analysis of dust emission complements 
dust extinction work because the dust emission method becomes most effective in the high 
column density inner regions of dense cores where dust extinction becomes large and difficult 
to penetrate. 

Much of the detailed information on starless core structure comes from observations of 
dust emission, using data from bolometer cameras (e.g. Ward- Thompson, Motte & Andre 
1999, Shirley et al. 2000, Visser, Richer & Chandler 2001). An important conclusion from 
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these studies is that starless cores appear to show flat density proflles in their inner regions, 
with extended envelopes that fall off rapidly in power law fashion. Recently, the predomi- 
nance of the flat central density gradients has been called into question as more sophisticated 
analysis including self-consistent temperature calculations indicate much smaller regions of 
flattening, or no flattening at all, in large part because the cores are cooler in their deep 
interiors than assumed previously (Evans et al. 2001, Zucconi et al. 2001). In addition, the 
central regions of even the nearest protostellar cores are generally comparable in size to the 
beamwidths of these telescopes and are poorly resolved. The density structure at smaller 
scales can be probed with interferometers. Strong constraints can be derived by analyzing 
interferometer data from dense cores directly in the visibility domain, though this approach 
has been rarely used (Keene & Masson 1990, Hogerheijde et al. 1999, Harvey et al. 2003a, 
Looney, Mundy & Welch 2003). 

L694-2 has been mapped at 850 //m with SCUBA on the JCMT by Visser (2000) and at 

1.2 mm with the IRAM 30m by Tafalla et al. (2003). The two maps have similar resolutions, 
14" for SCUBA and 11" for IRAM. The dust emission proflles suggest a steep outer density 
gradient {p ~ 2.7 from SCUBA, p = 2.5 from IRAM) consistent with that inferred from 
dust extinction (Harvey et al. 2003b). The usual isothermal analysis suggest a flattening 
of the density gradient within a radius of a few thousand AU. Visser (2000) flts a broken 
power law model that suggests an inner power law index of p = 0.8 within a break radius of 
8000 AU (~ 30"). Although no radio point source has been detected in the L694-2 core (e.g. 
Harvey et al. 2002), the measurement of the inner power law index suffers a potentially large 
systematic error from possible point-source contamination to the central beam (Shirley et 
al. 2002). In this paper, we present observations of dust continuum emission from L694-2 at 

1.3 mm obtained with the IRAM Plateau de Bure Interferometer (PdBI) and the Berkeley- 
Illinois-Maryland- Array (BIMA). The long baseline observations from the PdBI limit the 
contribution to the flux from any embedded point source, while the BIMA observations 
sample the dense core structure on size scales from ~ 10000 AU to ~ 1500 AU (resolution to 
6"). Together, these observations provide conclusive evidence of a flattening in the density 
profile within ~ 30" from the center of the core. 



2. Observations 
2.1. IRAM PdBI 

Continuum emission from L694-2 was observed at 1.3 mm (231.32 GHz) with the IRAM 
PdBI in the compact D configuration on 2001 November 20 (4 antennas) and 2002 April 03 
(6 antennas). Table 1 lists the observational parameters. A single pointing was used, with 
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center 19:41:04.44, 10:57:00.9 (J2000) chosen to coincide with the position of peak emission 

in the N2H+(l-0) spectral hne measured at BIMA (Jonathan WiUiams, private communica- 
tion). The half-power field of view for PdBI at this wavelength is 22" (5500 AU). The PdBI 
observations provide measurements at baseline lengths from the shadowing limit of 15 m to a 
maximum of over 100 m (resolution ~ 2'.'6). The absolute flux scale was set by observations of 
the standard source MWC 349, assumed to be 1.71 Jy. The estimated uncertainty in the flux 
scale is roughly 20%. Frequent observations of nearby calibrators J1751-I-096 and J1925-I-211 
were used to determine time-dependent complex gains. Continuum visibility records were 
formed for each 60 s integration of the digital correlator (3 x 160 MHz bandwidth, tuning 
double sideband). The correlator bandpass was measured with observations of the strong 
sources 3C 345 (2001 November 20) and 3C 273 (2002 April 03). The data were calibrated 
using the IRAM software package CLIC, and comprise a total of 7240 records (2040 records 
from 2001 November 20, 5200 from 2002 April 03). In addition to amplitude and phase, each 
record contains a variance measure, determined from the system temperature and antenna 
gains. 



2.2. BIMA 

Continuum emission from L694-2 was observed at 1.3 mm (231.32 CHz) with BIMA 
in the compact D configuration on 2001 September 22 and 29, 2002 June 08, and 2002 
September 09, and in the more extended C configuration on 2002 March 31. Table 1 lists the 
observational parameters. A single pointing of the nine antennas was used on each occasion. 
The pointing center for the BIMA observations is displaced from the PdBI pointing, by 
(5R.A.= -I'.'S, (5Dec.= 5'.'2. The half-power field of view for the BIMA antennas is 50" 
(12,500 AU). The BIMA observations provide measurements at basehne lengths from the 
shadowing limit of 6 m to a maximum of 70 m (resolution ~ 3'.'7). The bandpass and 
fiux calibration were determined through observations of Uranus, using the a priori antenna 
gains. The estimated uncertainty in the fiux scale is roughly 15%. Frequent observations 
of J1925+211 were used to determine time-dependent complex gains. Continuum visibility 
records were formed for each 23 s integration of the digital correlator (700 MHz bandwidth). 
The data were calibrated in the BIMA software package MIRIAD. The resulting ~ 2 x 10^ 
records were then averaged in time bins of 3.45 minutes (9 records). This choice of binning 
was made to reduce the number of visibilities to a more manageable size for analysis, without 
introducing significant phase error at the longer baselines. The resulting dataset contains 
22327 records. As with the PdBI dataset, a variance measure for each visibility measurement 
is also recorded. 
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3. Constructing Model Visibilities 

The visibility measurements are analyzed directly, without producing images that are 
limited by standard Fourier inversion and deconvolution techniques. This approach is compu- 
tationally intensive, but it allows a much more direct comparison with models than analyzing 
images. In particular, the results are not compromised by problems with the synthesized 
beam characteristics and dynamic range. 

The L694-2 visibilities are compared to theoretical models of protostellar envelope struc- 
ture by constructing synthetic visibilities, taking account of (1) the dust continuum radiative 
transfer, and (2) the specifics of the observations, including the exact {u, v) sampling and 
primary beam pointing and attenuation for the two telescopes. The models necessarily in- 
clude an assumed temperature distribution and (constant) specific mass opacity in addition 
to the model density field. We do not consider an exhaustive list of starless core models 
but instead fix attention on a few widely promoted density fields, including Bonnor-Ebert 
spheres, Phimmer-like models, broken power law descriptions, and isothermal cylinders. 
These models can all match the steep density gradient inferred for the outer regions of the 
core. 

The model datasets are constructed using the recipe described in Harvey et al. (2003a). 
In brief, a 512 x 512 model intensity image of resolution 0'.'5 pixel^^ is calculated using 
the full Planck function for the emissivity and integrating the radiative transfer equation 
through the model globule. Each model is normalized to a flux at 1.2 mm of 800 ± 80 mJy 
within a circular aperture (top-hat) of radius 30", calculated from the continuum map of 
L694-2 made by Tafalla et al. (2003) with MAMBO on the IRAM 30m (assuming a dust 
opacity spectral index of unity). We adopt an outer boundary of i?out — 0.15 pc in the 
models based on the extinction observations (Harvey et al. 2003b), although the results are 
not sensitive to this assumption. Observations are simulated by performing an EFT, and 
assuming a Gaussian form for the primary beams (PdBI EWHM 22", BIMA EWHM 50"). 
The exact {u, v) sampling is achieved by interpolating the real and imaginary parts of the 
resulting visibility grid. The center for each model is assumed to be at the pointing center 
of the PdBI observations, and in Section 5.2.1 we discuss the shght sensitivity of the results 
to this assumption. 

3.1. Model Selections 

As detailed above, the model visibilities are derived from the assumed (1) mass density 
distribution p(r), (2) dust temperature distribution Trf(r), and (3) specific mass opacity of 
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the dust Ki.3 mm (normalized to a constant flux within a 30" circular aperture). We consider 
the expected form of these quantities, and their uncertainties. 

3.1.1. Density 

The main goal is to constrain the density distribution of L694-2, given a realistic choice 
for the temperature distribution, and specific mass opacity. We consider four models of 
starless core density structure that can match the steep density gradient in the outer regions 
inferred from the near-infrared extinction. 

Bonnor-Ehert spheres. — These models are pressure-confined isothermal spheres, for 
which the solution remains finite at the origin (Ebert 1955, Bonnor 1956). They are solutions 
of a modified Lane-Emden equation (Chandrasekhar 1967): 

where ^ = [t/Rq) is the dimensionless-radius, Rq = a/ ^/ATiGp'c'is the (physical) scale-radius, 
and = ~^^{p/Pc) is a logarithmic density contrast, with pc the (finite) central density, 
and a is the efi^ective sound speed in the core (we adopt a = 0.20 km s~^ based on an 
assumed central temperature of 9 K, with turbulent component Oturb = 0.09 km s~^, as in 
the extinction study). Configurations with ^max > 6.5 are unstable to gravitational collapse. 
A highly supercritical dimensionless-radius of ^max = 25 ± 3 was found to reproduce the 
extinction observations of L694-2 for radii r > 30". 

Plummer-like models. — These are empirical models suggested by Whitworth & Ward- 
Thompson (2001) that capture the essential observed properties of starless cores with a 
minimum of free parameters. The model assumes that when a prestellar core becomes 
unstable against collapse at time t=0, it is static and approximates to a Plummer-like density 
profile (Plummer 1911), of the form: 

(2) 

The initial density is therefore uniform for r ^ i?0) and falls off as r^"^ for r ^ Rq. Whitworth 
& Ward-Thompson (2001) propose a fixed value of r/ = 4 in the model in order to reproduce 
the relative hfetimes and accretion rates for the Class and Class I phases. For L694-2, 
the density power law index in the outer regions inferred from near infrared extinction is 
p — 2.6±0.2 if the core is assumed to be unobscured, but has a steeper value of p = 3.7±0.3 if 
the core is embedded in a more extended uniform distribution of material, as is suggested by 



p(r, i = 0) = po 



Rn 



(Rl + r^y/^ 
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the flattening in the radial profile beyond r > 0.1 pc (Harvey et al. 2003b). A Plummer-like 
core with 77 = 4 embedded in a more extended cloud can therefore reproduce the extinction 
observations. The scale radius Rq is not constrained by the extinction measurements because 
no flattening is evident in the data (which do not penetrate the r < 30" inner region). 

Broken Power Law models. — These models represent a variation of the Plummer-like 
model, to allow for an inner density power law index that is non-zero. The density distribu- 
tion we adopt is of the form: 

^^""^ " {{r/Ro)(^P^ + {r/RQ)i'^ri)){i/2) (3) 

For r <^ Rq, the density falls off as r~^, and for r ^ Rq, the density falls off as r~'^. We 
choose this continuous prescription to prevent discontinuities in the flrst derivatives (which 
show up in the visibility profile) that would result from including a unphysical sharp break 
in the density distribution. We choose 77 = 2.6 to reproduce the observed extinction without 
introducing additional structure components. There is some degeneracy between Rq and 
p in the resulting visibility (or intensity) profile, since decreasing Rq or increasing p both 
produce a steeper emission profile or a fiatter visibility profile. The visibility amplitude for 
a typical starless core profile falls rapidly and is already very low at baselines that are long 
enough to sample the inner structure that cannot be probed with extinction. Even the high 
quahty visibihty dataset that we have obtained from BIMA and the PdBI therefore does 
not provide sufficient sensitivity to constrain both parameters simultaneously. We therefore 
choose a fixed value of Rq = 30" in order to quantify the effect of the inner power law index 
p. This choice is made on the basis that Rq = 30" is the largest value of Rq that is consistent 
with the extinction data, which will lead to the largest possible value of p allowed by the 
dataset (since a smaller value of Rq will necessarily require a shallower p to produce a given 
visibility amplitude at a given baseline). 

Cylinder models. — The isothermal cylinder model is a two dimensional analog of the 
Bonnor-Ebert sphere. The density is a function of the radial coordinate only (Ostriker 1964): 

where H — a/ y/AnGpc is the scale height that is equivalent to the scale radius Rq in the B-E 
analysis. The density is uniform near the axis of the cylinder but decays ever more rapidly 
with increasing radius, asymptoting to a power law of index p = 4 for r ^ H. The filament 
is supported radially by pressure gradients, but is unstable in the direction along its axis. 
The density distribution of the isothermal cylinder is a two dimensional case of the Plummer- 
like model, with \/8H = Rq, and with a physical basis for the normalization of the density 



-9- 



profile. A cylindrical model can therefore reproduce the observational properties of pre- 
stellar cores that provided the motivation for the Plummer-like models. In addition, because 
the isothermal cylinder and the B-E sphere both represent equilibria between self-gravity and 
gas pressure, the spherically averaged density profile of an isothermal cylinder can also mimic 
closely that of a B-E sphere, in particular a flat inner region with a steeply falling envelope 
(Boss & Hartmann 2001). The cylindrical model also provides a basis for interpreting the 
departures from spherical symmetry in the L694-2 core (Harvey et al. 2003b). A slightly 
tilted, embedded cyhnder with scale height H — 13.5"±1.5" reproduces the extinction profile 
for the inner 83" (0.1 pc) of the core. In the present study we consider isothermal cylinders 
viewed along the axis, since the subtle effect of a small tilt angle can not be constrained with 
the visibility dataset. 



3.1.2. Extended Cloud Structure 

Both Plummer-like and Isothermal Cylinder models can successfully describe the near- 
IR extinction observations, if the cores are embedded in an extended distribution of gas. 
This additional component to the density structure is suggested by the shape of the radial 
extinction profile, and must be accounted for in the fitting of these types of model. We include 
the extended cloud structure in the visibility analysis as follows. The extended structure 
is assumed to be smooth and is therefore resolved out by the interferometers. The only 
effect of the additional structure is to reduce the total flux that is attributable to the core in 
this context. The extinction profile of the L694-2 core asymptotes to a color-excess that is 
approximately 0.2 magnitudes higher than the background, and about a tenth (1/10) of the 
color excess at 30" radius (the edge of the flux normahzation aperture) . We approximate the 
intensity profile to be flat within 30" of the core, and assume that the temperature of the 
extended gas is equal to that in the core, so that the extended structure accounts for 10% 
of the flux normalization for each of these types of model. 



3.1.3. Temperature 

A detailed study of the expected dust temperature distribution in starless cores has 
been performed by Evans et al. (2001). They calculate the temperature distribution, Td(r), 
self-consistently using a ID radiative transport code, and assuming Ossenkopf & Henning 
(1994) opacities for grains that have grown by coagulation and accretion of thin ice mantles. 
They find that the Interstellar Radiation Field (ISRF) dominates the heating in the core, 
roughly a factor of 3 stronger than heating due to cosmic rays, even at the center of an 
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opaque core. The effect of the opacity law on Td{r) is small; the largest difference is at the 
center of the core, where opacities for coagulated grains that lack mantles cause a lower value 
of Td by ~ 0.5 K. 

Evans et al. (2001) present the dust temperature distribution Tfi{r) for a Bonnor-Ebert 
sphere with outer radius i?out = 0.17 pc (35000 AU) in which the gas is isothermal at 10 K 
(their Figure 3). The model has a central density of n = 1 x 10^ cm~^, and a dimensionless 
outer radius of ^max — 42, and is heated by an ISRF that combines the infrared behavior from 
COBE (Black 1994) and the ultraviolet behavior from Drainc (1978). The dust temperature 
varies from a minimum of ~ 8 K in the inner thousand AU of the core, to a maximum of 
~ 14 K in the outer regions. The model is similar in radius, but somewhat more centrally 
condensed than the B-E sphere that best fits the extinction in L694-2 (i?out — 0.17 pc, 
Ro = 3.4" in the computed model; R^ut = 0.15 ± 0.14 pc, Ro = 4.9" ± 0.6" for L694-2). In 
addition, the Evans et al. (2001) analysis neglects any shielding by an extended component 
of material that surrounds the core, as is seen around L694-2 in near-infrared extinction 
(Harvey et al. 2003b). The degree to which the interior of the L694-2 globule is shielded 
from the ISRF will differ somewhat from in the Evans et al. (2001) model. However, the 
Evans et al. (2001) model likely provides an accurate description of the actual temperature 
profile in L694-2 to within the uncertainties that stem from the inhomogeneities of the ISRF. 
We therefore adopt the Evans et al. (2001) temperature distribution for use in our modeling. 
We investigate the sensitivity to the assumed temperature distribution by also exploring 
models that assume a constant dust temperature of — 12 K. 

The Bonnor-Ebert and Cylinder models used in the model fitting are based on hy- 
drostatic equilibrium density configurations for a gas cloud that has an isothermal kinetic 
temperature T^. The Evans et al. (2001) temperature distribution is calculated using a 
density model that obeys the same assumption. At high densities efficient gas-dust coupling 
forces Tk to equal the dust temperature T^, while at low densities (n < 10^ cm"^) ^ 3^ 
(Takahashi, Silk, & Hollenbach 1983, Doty & Neufeld 1997). The Bonnor-Ebert model that 
best fits the L694-2 extinction has central density n ~ 3 x 10^ cm~^, with a center-to-edge 
density contrast of ~ 400. In the inner regions of L694-2 the kinetic temperature of the 
gas should decrease in a similar fashion to the dust temperature, but in the outer regions 
the two temperature distributions may depart. An entirely self-consistent approach to the 
problem would include a full calculation of the gas energetics, including dust coupling to 
calculate simultaneous density and (dust) temperature distributions. Such an approach has 
been followed by Galli, Walmsley & Gongalves (2002). However, to zeroth order, small 
departures from isothermality in the kinetic temperature have little effect on the density 
profile, since the pressure gradient of the gas is dominated by the gradient in the density 
and not the temperature. The effect of non-isothermality on the resulting emission profile is 
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therefore dominated by the dust temperature dependence of the Planck function. Since the 
uncertainties that surround the calculation of the dust temperature profile are substantial, 
the computational overhead of a more self-consistent approach is not justified at this time. 
Wc separate the problem into two parts, modeling density profiles assuming a constant gas 
kinetic temperature, but relaxing the isothermal assumption in order to model the resulting 
thermal emission from the core. 

A remaining issue deals with the appropriate temperature distribution to use for the 
cylindrical models, since these arc by definition not spherically symmetric. Since no self- 
consistent 3D dust radiative transfer results are available for this type of model, wc have 
made a very simple approximation, assuming the filament to have an aspect ratio of 2:1 and 
stretching the Evans et al. (2001) temperature distribution appropriately along the line-of- 
sight. 

3.1.4- Mass Opacity 

The mass opacity of dust grains in the millimeter region of the spectrum in protostellar 
envelopes is uncertain but generally assumed to follow a power law with frequency, Ki, oc u^. 
The power-law index varies depending on the dust properties, but tends to be bounded by 
a small range, roughly 1 to 2 (Ossenkopf & Henning 1994). The opacity of the dust affects 
the temperature distribution of the dust, as already described. For our purposes, only the 
spectral index is of importance, since we assume a constant opacity (independent of r), and 
the overall normalization is fixed by matching the MAMBO single dish observation. The 
spectral index affects the normalization of the models only weakly since the MAMBO fiux 
constraint is made at a wavelength very close to that of the PdBl and BIMA observations. 
We have adopted P = I; the extreme alternative of /3 = 2 leads to only a 4% change in fiux 
normalization (for Td = 12 K), which is small compared to the overall 10% uncertainty in 
the normalization itself. 



4. Method for Fitting Model Parameters and Evaluating Fit Quality 

The basic procedure is to maximize the probability distribution: 

P(M0del I data) = g-^^^-^^^-^'"^))'/"'^'' ^-{m-mo)V2al 
i 

where the Zi are the visibility data points with uncertainty (jj, f{xi;p,m) are the model 
data points, p a free parameter in the models, and m a model parameter about which 
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we have a constraint (namely that it is a Gaussian random variable with mean mo and 
standard deviation am)- Maximizing the probability distribution is equivalent to minimizing 
the logarithm of its inverse. Taking account the fact that the Zi are by nature complex 
visibilities, we want to minimize a modified x^: 



The sum/product can extend over any suitable subset of the visibility points. It is useful 
to be more explicit about the parameters p & m. The free parameter p is used to describe 
the shape of the model, e.g. the index of the inner part of the broken power-law density 
distribution. The parameter m allows us to include the observational uncertainties, e.g. 
the ~ 10% uncertainty in the normalization of the models derived from the Tafalla et al. 
(2003) map, and the ~ 15% uncertainty in the fiux calibration for the BIMA data. This is 
achieved by allowing the BIMA model visibilities to be scaled by a constrained parameter 
m, assumed to be a Gaussian random variable with mean rriQ = 1.0 and standard deviation 
(Jm = 20%. The longer baselines of PdBl do not detect any signal from the L694-2 core. 
These measurements are used to constrain a limit on the point source contamination and 
are not included in the model fitting. 

For a given model, we evaluate the best fit parameters by minimizing the modified 
distribution. Uncertainties in the parameter vahics arc analyzed using the Monte Carlo 
technique known as the bootstrap (Press et al. 1992). In brief, the dataset is resampled n times 
(typically n ~ 200), each time the fitting process is repeated and the best-fit parameters 
recorded, until the distribution of best-fit parameters is well sampled. The width of the 
distribution provides an estimate of the uncertainty in the parameters that best fit the 
original dataset. 

This numerical approach is made necessary by the non-linear nature of the fitting pa- 
rameters, and is especially useful for this analysis because small variations in the value of the 
modified ^ defined above may represent surprisingly large variations in fit quality. The 
nature of the visibility dataset, comprising a very large number of very low signal-to-noise 
measurements, means that even a model with zero signal will on average reproduce each 
visibility measurement to within its uncertainty. The resulting shallowness of the wells 
causes two models with slightly differing parameters to seem almost equally good despite 
the fact that there is ample signal to distinguish them. The numerical approach essentially 
bypasses this comphcation. The issue could also be circumvented by binning the visibilities 
(e.g. radially in {u, v) distance), to increase the signal-to-noise, and then performing a fit 
to the binned values. We avoid this solution because the visibility profiles for the various 
models are falling steeply at the baselines of interest, and binning the data inevitably intro- 
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duces a bias in the fit due to contracting the range of basehne vectors to a single length at 
the center of the bin. We use binning only as a graphical tool in order to demonstrate the 
fit quality of a particular model. 

5. Results and Analysis 

5.1. Limits on Point Source Flux 

The PdBl datasct covers baselines from 12 to 80 kA, with an overall rms noise of 0.9 mJy. 
The visibility data are consistent with the noise. The long baselines covered by the PdBI 
data therefore provide an important constraint on the maximum compact component of the 
flux from the L694-2 core, F^. < 2.7 mJy (3 a). For optically thin dust, this limiting flux 
corresponds to an implied (disk) mass limit of M < 5 x 10~^ Mq (60 K/Tdisk) for an opacity 
K'l.a mm = 0.02 cm^ g~^. This hmit is roughly an order of magnitude lower than the 0.002- 
0.3 M0 range of disk masses observed around T-Tauri stars by Beckwith et al. (1990), and 
further demonstrates the "starless" nature of the L694-2 core. 

The flux limit allows an estimate of the maximum bolometric luminosity for any compact 
component embedded in the L694-2 core. For a simple estimate, we assume a dust opacity 
spectral index of unity, and model the compact component as a graybody of the form: 

= 5,((Td,st)) (1 - exp[-r,]) , (7) 

where -By((Tdust)) denotes the Planck function at frequency u for a mean dust temperature 
(Tdust), is the dust optical depth, and fls the solid angle subtended by the source (e.g. 
Beckwith et al. 1990). Since the envelope is entirely optically thin at 1.3 mm, the graybody 
must have a flux that is < 2.7 mJy at this wavelength. For a given mean dust temperature, 
this constraint fixes the mass of the compact component. The only remaining parameter 
is the radius R of the component, which essentially identifies the wavelength at which the 
emission becomes optically thick. 

For a compact component to remain undetected in our PdBl observations implies that 
its bolometric luminosity is: 

Lboi < 0.07 Lo (i?/100 AU)°-i ((rdisk)/50 K)^-^ (8) 

Note that these power law dependencies are an approximation and are not accurate for 
changes in the parameters by factors of more than ~ 2. The dependences on temperature and 
radius can be understood in the following way. The dominant contribution to the luminosity 
comes from the R-J region of the Planck function, in which the emission is optically thin 
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(for this typical size scale and temperature, the compact component becomes optically thick 
slightly to the Wien side of the peak in the Planck function). If the emission from beyond 
the peak in the Planck function is entirely negligible, then the luminsosity is proportional to: 
L oc /(I'™'''' Kyv'^du. Note that the flux constraint at 1.3 mm causes there to be no temperature 
dependence in the integrand. Using the Wien law for z/max oc (T), and Ky oc z/, then gives 
Lboi oc (T)"^. The contribution of the optically thick emission causes the actual temperature 
dependence to depart from this relation, and introduces a weak dependence on the radius of 
the compact source. 

The limit the PdBl observations place on the maximum bolomctric luminosity of any 
embedded compact object is substantially lower than the IRAS limit for L694-2 of ~ 0.3 L©. 
Reprocessing of the warm compact emission to long wavelengths by the envelope essentially 
allows a large bolometric luminosity to remain hidden below the IRAS detection limit. 

5.2. Density Structure 

The PdBI data constrains the maximum compact component of dust emission, which 
allows us to proceed to fit for the envelope structure using the BIMA data at shorter base- 
hnes. The envelope visibility profiles fall steeply with increasing basehne. (Note that a 
fiatter emission profile leads to a steeper visibility falloff.) At the longer BIMA basehnes, 
the visibility amphtudes for the best fit models falls well below the signal-to-noise of the 
entire BIMA dataset. This makes restricting the visibility fitting to the shorter baselines 
advisable, since extending the fitting range to longer baseline adds noise but no signal. Hav- 
ing experimented with various upper limits on the baseline length, we found that an upper 
limit of 10 kA was optimal. Increasing the upper limit to 15 kA, or 20 kA did not change the 
results but produced higher values of the reduced x^. The visibility dataset to which the 
models were fit contains 6424 visibility measurements (each a 3.45 minute integration) and 
covers baselines from 5 to 10 kA. The rms noise in this reduced dataset is 3.3 mJy. Table 2 
summarizes the results from the fits to the four types of model described in Section 3.1. 
Below we describe these results. 

Bonnor-Ehert Sphere. — the model that best fits the visibility data has dimensionless 
outer radius ^max = 18 1| (Fit I). At the 2 a level this is consistent with the results of the 
extinction study, where a B-E sphere with ^max = 25 ± 3 reproduced the observed color 
excess. 

Plummer-like model. — the embedded model that best fits the visibility data has scale 
radius i?o = 26" or 6500175°° AU (Fit II). This turn-over radius is consistent with the 
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Visser (2000) analysis of SCUBA data, which fit a broken power law with a transition to a 
shallower index within a break radius of ~ 32" (8000 AU). 

Broken Power law. — the model that best fits the visibility data has a power law index 
in the inner region (r < 30") of p = O.Qloiie (^^^ m)- The turn-over radius used in this 
model is similar to that from the Visser (2000) analysis, and the fitted inner power law index 
is consistent with their result of = 0.8 (no uncertainties are given), despite their use of an 
isothermal dust temperature distribution. Note that the use of a smaller turn-over distance 
results in a lower fitted value of p. The above result should therefore be construed as a 
"mean" index in the r < 30" (7500 AU) region, or as an upper limit on the power law index 
in the innermost regions. 

Cylinder. — the end-on cylindrical model that best fits the visibility data has scale height 
H = 12" +5;° ", or 3000+^^JJ AU (Fit IV). This is consistent with the results of the extinction 
study, where a slightly tilted cylinder with H = 13.5" ± 1.5" best matched the observed color 
excess. 

For each of the best fit models, the fitted value of the scaling parameter is m ~ 0.9±0.05. 
That the preferred value is similar for each model and not unity suggests that the data prefer 
a particular slope of visibility profile that when extrapolated back to zero baseline has 10% 
lower normalization than the MAMBO flux. This difference is well within the combined 
uncertainties of the MAMBO flux measurement and the BIMA calibration accuracy (total 
~ 20%). 

Figure 1 shows plots of visibility amplitude vs. {u, v) distance for the various best fit 
models, and models that differ by ±2 a in the fitting parameters. The profiles of the best fit 
models are remarkably similar, and we return to this point in Section 5.3. The mean signal in 
the BIMA dataset used to fit the models is shown on the plots (filled circle) with ±2 a error 
bar. The binning introduces an uncertainty (and/or bias) in the appropriate {u,v) distance 
for the data point, that is not present in the model fitting procedure. For each type of model, 
we have adopted the {u, v) distance that provides a match with the best-fit model. This is 
motivated by purely graphical reasons, but we note that the mean baseline length weighted 
by the signal to noise of the best fit model (for constant noise on each baseline) provides 
an essentially indistinguishable {u, v) coordinate. However, the greater spread of the fitted 
models compared to the binned data point demonstrates the effect of this "smearing" of the 
signal in {u, v) distance. 

The plots in Figure 1 also show ±2 a confidence intervals of an azimuthally averaged 
synthetic visibihty profile derived from the Tafalla et al. (2003) MAMBO map of L694-2. 
The visibility profiles are computed assuming the MAMBO map to be the true intensity 
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distribution convolved with a Gaussian beam of 11" FWHM. No assessment has been made 
for systematic errors in the profile due to chopping and subsequent reconstruction of the 
MAMBO observations. The synthetic confidence intervals are extended until the the signal 
to noise drops below 2 (i.e. the intervals become consistent with zero signal). 

5.2.1. Systematic Uncertainties 

Table 3 lists the main sources of systematic uncertainty along with the level of un- 
certainty they produce in the fitted density parameters (.^max, ^o, P, & H). The various 
systematic errors are discussed below, in rough order of importance: 

Dust temperature distribution. — This is the largest source of systematic uncertainty. 
The temperature distribution that has been assumed in the analysis varies from ~ 8 K 
at the center to ~ 14 K in the outer regions, and it is motivated by physical argument. 
Nevertheless, the profile is probably uncertain at the level of ±2 K. To study the effect of 
this uncertainty on the inferred density structure, we have repeated the analysis for each 
model, replacing the physical temperature distribution with an isothermal distribution at 
12 K. The removal of the central minimum in the temperature profile results in a given 
density model producing a steeper emission profile and a shallower visibility profile. The 
density structure inferred is therefore less centrally concentrated than with the physical 
temperature profile. The effect is of a similar magnitude for all of the models: the best-fit 
parameters are changed by close to 2 cr. 

Unresolved compact component. — The PdBI data limit the maximum contribution to 
the emission from an unresolved compact component: F(. < 2.7 mJy (3 a). Including a 
dimmer compact component in the model fitting essentially reduces the visibility amplitude 
that is attributable to the envelope and results in fitted density structures that are more 
shallow. For a point source flux of Fc = 2 mJy, the density structure parameters are changed 
by roughly 1 a. 

Central position. — The flatness of the inner density profile and the lack of a detected 
compact component make the appropriate central position for the L694-2 core uncertain. 
However, the baselines of the BIMA visibilities from which the density structure is inferred 
are short and correspond to a fringe size of ~ 20". This means that the results are not 
very sensitive to small changes in the adopted central position (the pointing center of the 
PdBI observations). The signal in the combined visibilities peaks at 10 mJy within 1" of the 
adopted central position — the peak of N2H"'"(l-0) emission. A shift in the central position 
of 5", i.e. the shift between the fitted center from the extinction study and central position 
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assumed here, leads to a ~ 2 mJy reduction in the binned signal and a 1 cr change in the 
density structure parameters, corresponding to shallower density profiles. 

Extended structure. — A leading source of systematic uncertainty in the L694-2 extinc- 
tion study was the effect of the contribution of the extended structure in which the core 
appears to be embedded (Harvey et al. 2003b). Including this additional component to the 
models lead to inferred core density profiles that were significantly more steep, consistent 
with Plummer-likc and isothermal cylinder models. In the present study, this uncertainty 
has much smaller effect on inferred structure. Since the extended structure must be largely 
resolved out by the interferometers, neglecting the extended structure entirely leads to only 
a 10% increase in the amplitudes of the visibility profiles of the model cores. Repeating the 
fitting analysis for the two types of model in this context results in inferred density profiles 
that are ~ 0.5 u more shallow. 

Outer boundary. — The outer boundary of the L694-2 core has little effect in the model- 
ing due to the hmited field of view of BIMA and PdBI. For BIMA, the intensity distribution 
at the assumed outer edge of the core is suppressed by a factor ~ 0.05 due to the antenna 
pattern. Moreover, the intensity at the edge is already lower than the intensity at the center 
by a factor of ~ 500 (Fit I). Modifying the outer boundary therefore has no effect on the 
shape of the inferred density profile. However, for the Bonnor-Ebcrt sphere (Fit I), changing 
the outer radius docs modify the density structure parameter ^max because this parameter 
describes the dimcnsionlcss profile, not the physical profile. For this fit, the alternative 
structure parameter, the scale radius Tq = -Rout/Cmax is unchanged. 

Dust opacity spectral index. — The dust opacity adds systematic error via the use of 
the spectral index to transform the fiux normalization from MAMBO to the frequency at 
which our observations were made. To some extent, the appropriate MAMBO frequency 
itself depends on the spectral index, because of the broad bandwidth used in the single 
dish bolometer observations. The spectral index of the dust opacity should lie in the range 
P —1-2. We have used P — I throughout; an alternative extreme, /? = 2, leads to only a 4% 
reduction in the flux normalization, and to inferred density structure that is more shallow, 
but to a degree that is negligible in comparison to the other sources of error (both systematic 
and random). 

5.3. Discussion: A Physical Density Distribution for L694— 2 

The fitted physical models demonstrate clearly that the steep density gradient observed 
in the r > 30" region with extinction (Harvey et al. 2003b) does not continue in the inner 
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region that could not be probed in that study. This is consistent with the modehng of single 
dish dust emission observations by Visser (2000) and Tafalla et al. (2003). 

The visibility profiles of the best fit models shown in Figure 1 are remarkably similar to 
each other, only the power law model differing enough to be distinguishable by eye, and then 
only at long baselines. This provides an interesting demonstration of the degeneracy of the 
various models for starless core density structure. Figure 2 presents a plot of the best-fit den- 
sity models. The normalization of the density assumes an opacity of Ki. 3^111 = 0.02 cm^ g^^, 
which is uncertain by a factor of ~ 5 or more (Ossenkopf & Henning 1994). At large radii the 
Bonnor-Ebert and Power-law models differ from the Plummer-like and cylinder models due 
to the fact that the latter models are assumed to be embedded in an extended uniform dis- 
tribution of gas. The Plummer-hke and Bonnor-Ebert models are almost identical, and have 
central densities that differ by only 5% (n(H2) = 1.4 x 10^ cm~^ with uncertainty ~ 50%). 
The density profile of the end-on cylinder departs significantly from the the Bonnor-Ebert 
and Plummer-like profiles. The profile is more shallow with much lower central density 
(n(H2) = 5.1 X 10"^ cm~^ with uncertainty ~ 30%). This occurs despite the identical form 
of the expressions for the cylinder and Plummer-like models (with Rq = ^/8H), and the fact 
that the cylinder and Bonnor-Ebert models both represent a balance between self-gravity 
and thermal pressure. The difference derives from the lower dimensionality of the cylindrical 
model; a given line-of-sight corresponds to a constant "radius", and hence there is no ra- 
dial integration that for the spherically symmetric models makes the column density profiles 
more shallow than the density profiles. In addition, the extension of the cylinder along the 
line-of-sight means that the densities are correspondingly lower at a given radius. The asym- 
metry of the L694-2 core viewed in extinction provides a basis for preferring the cylindrical 
model over the others. A tilted cylinder with H = 13.5, Lsin0 = 0.14 pc {L ~ 0.2-0.5, and 
~ 20-45°), and central density (n(H2) = 4 x 10'^ cm^^ embedded in a uniform distribution 
of gas with column density N{H + H2) ~ 6 x 10^^ (L/0.5 pc) cm~^, successfully reproduces 
the dust emission visibility profile, as well as the profile and asymmetry of the dust extinction 
map. The instability of the cylindrical model along its axis is also consistent with the inward 
motions in L694-2 inferred from molecular spectral fines (Lee, Myers & Tafalla 2001). 

As already noted, the inner power law index of the fitted broken power law model 
(p = 0.9 to.ll: Fit II) iiiay be construed both as a "mean" index in the r < 30" region, and as 
an upper limit on the index in the innermost regions (due to the trade-off between turn-over 
radius and power law index in the inferred visibility profile) . The former is illustrated by the 
fact that the radially averaged "mean" effective power law index for the best-fit Plummer- 
like model is p = 0.86 over the range in radius < r < Rq (with Rq = 26"). The profile 
is consistent with the Visser (2000) study who found that a similar index {p = 0.8) and 
turn-over radius {Rq — 32") reproduced the SCUBA map. The index in the inner region is 
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much less than that observed in the outer envelope and shows that the density distribution 
in the inner region of the L694-2 core has clearly not relaxed fully, along the lines assumed 
in the inside-out collapse of Shu (1977), despite the presence of extended inward motion in 
the gas (Lee et al. 2001). These inward motions may reflect the contraction of a prolate core 
along its major axis (Harvey et al. 2003b). 

6. Summciry 

We present a study of the density distribution of the candidate contracting starless 
core L694-2 using high resolution 1.3 mm dust continuum observations from BIMA and the 
IRAM PdBI. In summary: 

1. The PdBI visibility data span baselines from 12 to 80 kA (spatial scales from ~ 500 to 
~ 5000 AU) and do not detect any significant emission from the core. This provides 
a stringent constraint on the maximum point source flux: Fc < 2.7 mJy (3 a). This 
flux hmit corresponds to a very small disk mass, M < 5 x 10~^ M© (60 K/Tjisk), and 
bolsters the "starless" interpretation of the L694-2 core. 

2. The BIMA visibility data in the baseline range 5 to 10 kA (spatial scales ~ lO'^-lO'^ AU) 
detect emission from the core and are used to constrain models for the density structure. 
That no signal is evident on baselines beyond 10 k A confirms the fiat emission profile 
in the inner regions, and a significant change in behavior from the steep density profile 
for r > 8000 AU inferred from near infrared extinction of background stars (Harvey et 
al. 2003b). We fit four types of model for starless core density structure that can match 
the steep density gradient in the outer regions inferred from extinction. We adopt a 
temperature distribution that decreases in the inner regions due to shielding from the 
ISRF (Evans et al. 2001). The best fit Bonnor-Ebert sphere has dimensionless outer 
radius ,^max = 18 ^'^4 (1 a). The best fit Plummer-like model has turn-over radius Rq = 
26 "tl " (6500^^^^° AU). The best fit broken power law model has index p = 0.9 toU 
(this model uses the maximum allowable turn-over radius of Rq — 30" (7500 AU) and 
provides an upper limit on the index at the center of the core). The best fit end-on 
isothermal cylinder has scale height H = 12 "tH " (SOOO^J^g AU). 

3. We consider the effects of various sources of systematic uncertainty on the derived den- 
sity structure. The largest uncertainty comes from the dust temperature distribution; 
assuming an isothermal core modifies the best fit parameters by ~ 2 cr, and makes the 
inferred density profiles more shallow. A possible weak compact component to the 
emission (F ~ 2 mJy) and the appropriate central position {59 < 4") are the main 
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remaining sources of uncertainty. Both of these uncertainties lead to ~ 1 cr changes in 
the inferred model parameters, again to shallower density profiles. 

4. The density profiles of the various best-fit models are nearly indistinguishable at base- 
lines shorter than ~ 10 kA. There is a strong degeneracy between the models for 
starless core density structure. The two physical models (B-E sphere and isothermal 
cylinder) represent a balance between self-gravity and thermal pressure, and therefore 
they produce similar profiles. The asymmetry of the L694-2 core evident in the extinc- 
tion map may provide a basis for favoring the cylindrical model. A cylinder slightly 
tilted to the line-of-sight with H = 13.5", Lsin0 = 0.14 pc (L ~ 0.2-0.5, and ~ 20- 
45°) and central density {n{}l2) = 4 x 10^ cm~^ embedded in a uniform distribution of 
gas with column density N{H + H2) ~ 6 x 10^^ (L/0.5 pc) cm~^ reproduces both the 
dust emission and extinction measurements. 
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Table 1. Summary of instrumental parameters 



Parameter 


PdBI Value 


BIMA Value 


Observation dates 


2001 Nov. 20; 


2001 Sep. 22; Sep. 29; 




2002 Apr. 03 


2002 Mar. 31 Jun. 08; Sep. 09 


Configuration 


D (four ants.); D (six ants.) 


D; D; 0; D; D (nine ants.) 


Pointing center (J2000) 


19''41™04'.'44, +10°57'00'.'9 


19''41'"04'.'32, +10°57'06'.'l 


Observing frequency 


231.3 GHz 


231.3 GHz 


Phase calibrators 


J1751+096, J1925+211 


J1925+211 


Bandpass calibrator 


30 345; 30 273 


Uranus 


Flux calibrator 


MWO 349 


Uranus 


Primary beam FWHM 


22" 


50" 


Bandwidth 


480 MHz 


700 MHz 


RMS 


0.9 mJy 


1.6 mJy 





Table 2. 


Summary of the Density Model Fits 




Fit 


Density Model 


Fitted Model Parameter 


Oalibration scaling 




I 

n 

HI 
IV 


Bonnor-Ebert 
Plummer-like 
Broken Power Law 
Oylinder 


^max = 18 ^4 

Ro = 26"tt" 

p - u.y -0.16 
H = 12" Hi " 


M = 0.90 ±0.04 
M = 0.90 ±0.05 
M = 0.92 ±0.04 
M = 0.90 ± 0.05 


1.2520 
1.2521 
1.2519 
1.2520 



-24- 



Table 3. Summary of systematic uncertainties in fitted parameters 



Model Assumption 


Variation 


Resulting Systematic Error (a) ^ 


Temperature distribution 


Td = Tphys(r) ^ 12 K 


~ 2 


Central point source flux 


F ~ 2 mJy 


~ 1 


Central position 


se < 4" 


~ 1 


Neglect extended structure ^ 


Fo{r < 30") = 80 mJy — > 


~ 0.5 


Outer boundary of L694-2 ^ 


SRout/Rout < 50% 


« 1 


Dust opacity spectral index 


(3=1 — >2 


«: 1 



^Systematic error in the fitted model parameter given in units of a the random error in the 
parameter listed in Table 2. Note that each quoted systematic error corresponds to a shallowing 
of the inferred density profile. 

^Applies to Plummer-like and Cylinder models only 

^Does not apply to B-E sphere; for the B-E model, increasing the outer radius increases ^max 
in proportion to .Rout, keeping constant scale radius ro 
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Fig. 1. — Visibility amplitude vs. (u.v) distance for L694-2 at 1.3 mm, for the best-fitting 
models (dashed lines) and models that deviate by ±2 a in the fitting parameters (dotted 
fines). The binned BIMA visibifity ampfitude in the baseline range 5 to 10 kA is also shown 
(circle), with a ±2 a error bar. The PdBI limit on a point-like component (> 12 kA) is shown 
as a horizontal (dash-dot) fine. The sofid lines represent ±2 a intervals of the azimuthally 
averaged synthetic visibility profile derived from the MAMBO map. See text for further 
explanation. 
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Fig. 2. — Number density of molecular hydrogen vs. radius for the various best fit models of 
L694-2. The number density calculation assumes a dust opacity of Ki.smm = 0.02 cm^ g~^, 
a mean molecular weight of fi = 2.29, and a Hydrogen mass fraction Xji = 0.73. The 
best-fit Bonnor-Ebert (dotted) and Plummer-like (dashed) profiles are nearly identical. The 
cylindrical model (solid) has density a factor ~ 2 lower due to the extension of the cylinder 
along the line-of-sight. 



